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Abstract. We discuss the regular magnetic field in the 
galaxy M51 on the basis of a new interpretation of po- 
larization angles observed at the wavelengths AA2.8, 6.2, 
18.0 and 20.5 cm. We found a magneto-ionic halo in M51 
with a radial extent of about 10 kpc. The regular magnetic 
fields in the disk and the halo have different structures. 
The regular magnetic field in the halo is axisymmetric 
and horizontal. Its field lines are spirals pointing inwards 
and generally opposite to those in the disk. 

The azimuthal structure of the magnetic field in the 
disk is neither axisymmetric nor bisymmetric but can be 
represented by a superposition of these two basic harmon- 
ics with about equal weights. Magnetic lines of the regu- 

■ lar field in the disk are spirals generally directed outwards. 
They are well aligned with the optical spiral arms as traced 

■ by the dust lanes. 

The regular magnetic field strength averaged in 3 kpc 
wide rings is about 5-10 fiG in the disk and 3 /J.G in 
the radial range 3-6 kpc in the halo. These values are in 
good agreement with independent estimates from the syn- 
chrotron emission. The general features of the magnetic 
patterns revealed in the disk and the halo seem to be in 
agreement with predictions of dynamo theory. 

Our results are also based on an analysis of the 
magneto- ionic medium in M51. For a range of radii in 
the galactic disk, we estimated the scale heights of syn- 
chrotron and thermal disks, electron densities and filling 
factors as well as the amount of depolarization. 

Key words: spiral galaxies - interstellar magnetic field 
- radio continuum - polarized radio emission - Faraday 
rotation - M51 - depolarization 
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1. Introduction 

M51 (NGC 5194) is the first external spiral galaxy 
from which linearly polarized radio emission was detected 
(at A20cm by Mathewson et al. 1972, and at AA6.0 and 
21.2cm by Segalovitz et al. 1976), and for which the global 
structure of the regular magnetic field highlighted by these 
observations was investigated (Tosa & Fujimoto 1978). It 
is also the first galaxy for which the discovery of the so- 
called bisymmetric magnetic structure was claimed (Tosa 
& Fujimoto 1978), a configuration of magnetic field which 
was later suspected for some other nearby spiral galaxies 
(see Krause 1990, Beck 1993 and Beck et al. 1996 for a 
review) . The explanation of the origin of the bisymmetric 
magnetic fields has become one of the major challenges to 
the theory of galactic magnetic fields of the last decade. 

Later observations of M51 with better sensitivity 
and resolution (Neininger 1992a; Horellou et al. 1992; 
Neininger et al. 1993a; Neininger & Horellou 1996), and 
also a more careful analysis of the observational data, 
revealed a magnetic pattern which is considerably more 
complicated than a simple bisymmetric structure. This is 
true also for some other galaxies, e.g. M83 (Neininger et 
al. 1991, 1993b). 

One of the goals of the present paper is to introduce 
and test a general way to represent the complicated mag- 
netic patterns observed in galaxies in terms of a tractably 
small number of parameters. Hopefully, this will facilitate 
a fruitful confrontation of theory with observations. Such 
a parametrization can be conveniently performed in terms 
of the Fourier expansion of the magnetic field in azimuthal 
angle. The lowest Fourier harmonic corresponds to the 
axisymmetric magnetic field, the next higher one to the 
bisymmetric mode, etc. However, we emphasize that the 
Fourier harmonics thus derived are not necessarily con- 
nected with the dynamo modes and their physical mean- 
ing should be established using models of the magnetic 
field generation and evolution in a given galaxy. 
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Fig. 1. Maps of the E-vectors rotated by 90° observed at AA2.8cm (a), 6.2cm (b), 18.0cm (c) and 20.5 cm (d). The length 
of the vectors is proportional to the observed polarized intensity. They are shown superimposed onto an optical picture (Lick 
Observatory) of M51. All maps were smoothed to an angular resolution of 75" 



We attempted to build a coherent, self-consistent pic- 
ture of the global magnetic structure based not only on the 
Faraday rotation analysis, but also on other available in- 
formation coming from, e.g., intrinsic polarization angles, 
depolarization data, total synchrotron emission, thermal 
radio emission, the morphology of the galaxy, etc. In order 
to interpret a polarization pattern in a galaxy, one should 



know certain parameters of the interstellar medium such 
as scale heights of the thermal and synchrotron disks, elec- 
tron volume density and filling factor. These are discussed 
in Sect. 3. 

On pursuing our goals we used recent multifrcquency 
observations of M51, which allowed us to distinguish two 
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magneto-ionic layers along the line of sight (the disk and 
halo) with significantly different magnetic fields. 

We adopted the following parameters of M51: centre 
coordinates a 50 = 13 h 27 m 46*327, 5 50 = +47 2 7'lf/.'25 
(Ford ct al. 1985), a position angle of the major axis of 
— 10° measured counterclockwise from north, an inclina- 
tion angle i = —20° (i = 0° is face-on - see also Ap- 
pendix A in the electronic version) (Tully 1974), and a 
distance to M51 of 9.7 Mpc (Sandage & Tammann 1974). 

In this text equations, figures and tables are numbered 
as in the electronic version, whereas the sections are num- 
bered consecutively. 



2. The observational database 

The observations of the spiral galaxy M51 which we 
discuss below were obtained at the wavelengths 2.8 cm 
(Neininger 1992a), 6.2cm (Neininger et al. 1993a), 18.0cm 
and 20.5cm (Horcllou et al. 1992). In contrast to earlier 
discussions of these observations, we analyze here the ob- 
servations for all four wavelengths simultaneously. 

The 2.8 cm data are single-dish measurements ob- 
tained with the 100-m Effclsberg radio telescope. The 
other data sets were obtained with the VLA in its D- 
array. This imposes some restrictions for the use of the 
6.2 cm data since at this wavelength the diameter of the 
primary beam of the VLA is 9 arcmin, which corresponds 
to a diameter of 25 kpc in the plane of M51. Therefore, we 
considered the measurements at 6.2 cm to be reliable up 
to a radius of 9 kpc chosen to be somewhat smaller than 
the radius of the primary beam, and we did not use them 
at larger radii. 

At all wavelengths the data were smoothed to a final 
resolution of 75" corresponding to 3.5 x 3.8 kpc in the plane 
of M51. Fig. 1 shows the observed E- vectors rotated by 
90°, superimposed onto an optical picture of M51. 

Following a usual procedure in the studies of regular 
magnetic fields in external galaxies, the galaxy was di- 
vided into several rings and we considered the values of 
polarization angle averaged in sectors in each ring, tp ni , 
and the corresponding uncertainties, a n i, where the sub- 
script n refers to the wavelength and i to the sector. Here 
we chose the rings between the galactocentric distances 
r = 3, 6, 9, 12 and 15 kpc; within each ring sectors of an 
opening angle of 20° were used. The azimuthal angle 9 was 
measured counterclockwise from the northern major axis. 
Throughout the paper we specify sectors by their median 
value of 9. A detailed discussion of the data averaged in 
sectors and their uncertainties is given in Sect. 2.2 of the 
electronic version. 



3. The magneto-ionic medium in M51 

3.1. The nonthermal disk in M51 

The exponential scale heights of the synchrotron emission 
at A2.8cm and A20.5cm given in Table 1, h 2 .s and /120.5 
respectively, were estimated from those in the Milky Way 
by scaling the latter values obtained at 408 MHz (Beuer- 
mann et al. 1985) with frequency as i/~ - 25 , as observed 
for NGC 891 (Hummel et al. 1991a) and M31 (Berkhuijsen 
et al. 1991). As M51 and the Milky Way are galaxies of 
a similar type (Sc and Sbc, respectively) and have about 
the same linear dimensions, the scale height in the So- 
lar neighborhood (r Q = 8.5 kpc) was assumed to apply at 
r = 9 kpc in M51. 

The total field strength B (including the regular and 
turbulent components) can be evaluated from the non- 
thermal radio intensity per beam area at A2.8cm, J 2 .8, 
using the standard assumption of equipartition of energy 
density between the magnetic field and the cosmic rays 
(see, e.g., Krause et al. 1984). The nonthermal emission 
was obtained from the total emission by subtracting the 
thermal component derived by Klein et al. (1984). 

Using the total nonthermal intensity one first estimates 
the strength of the transverse magnetic field (i.e. the pro- 
jection of B on the plane of the sky). From the polarized 
intensity, or the degree of polarization, one obtains the 
strength of the transverse regular magnetic field _B rog± , 
from which one gets B reg by deprojection assuming that 
the field lies in the galaxy's plane. The turbulent magnetic 
field strength i?tur is found from its transverse component 
by multiplication by -\/3 assuming statistical isotropy. As 
Faraday effects are negligible at A2.8 cm, the degree of po- 
larization at this wavelength, P2.87 yields the best estimate 
for the strength of the regular magnetic field. Therefore 
we evaluated B, B reg and -B tU r from the A2.8cm data as 
given in Table 1. 

We stress that the estimates of B reg and B tm from the 
synchrotron emission are used only to assess the role of 
Faraday depolarization effects in Sect. 3.3. The analysis 
of the polarization pattern performed in Sect. 4 yields in- 
dependent estimates of B reg which are in agreement with 
those given in Table 1. 

3.2. The thermal disk in M51 

Klein ct al. (1984) derived the radial dependence of the 
thermal radio emission in M51 at A2.8cm. For each ring 
the average thermal flux density 6*2.8 per beam area is 
given in Table 2. To calculate the average electron den- 
sity we used the formulae for the thermal flux density (in 
terms of emission measure) of Hn regions of Mezger & 
Henderson (1967) in the form given by Israel et al. (1973) 
for an unresolved source. In our case we have 

K) = Cfl/2 ( s 2 . s \ 1/2 /2^\ - 1/2 ^ 

cm -3 \ mJy/beam J \ pc J 
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Table 1. The synchrotron disk in M51 
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Table 2. Thermal disk and halo in M51 
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Note: the negative error in £2.8 was cut off to prevent (n e ) < 0. 



where (n c ) is the average electron density in the ther- 
mal ionized gas layer of exponential scale height hth , f is 
the volume filling factor of the electron density defined by 
(n e ) 2 = f(nl), and C is a certain constant depending on 
distance, resolution and electron temperature. For a dis- 
tance of 9.7 Mpc, a resolution of 75", and a temperature 
of 10 4 K we have C = 7.4. 

Using the values of / and hth derived below we found 
the average electron densities presented in Table 2. 

The scale height hth is not known for M51. Again we 
can make an estimate using the values known for the Milky 
Way. Observations indicate that hth varies with radius. 

For the inner Galactic region we used the observation 
of Reich & Reich (1988) that the full halfwidth of the 
thermal emission at r < 7 kpc is 2? 0-2? 5, which yielded a 
scale height of about 400 pc at a mean radius of r — 4 kpc. 
This scale height was adopted for the ring 3-6 kpc in M51. 

For the solar neighbourhood Reynolds (1991a) sug- 
gested a two-component model of the thermal disk. For 
our present purpose it is sufficient to consider a one- 
component model of the thermal disk described by a sin- 
gle exponential fitted to the sum of the thin and the thick 
disk (see Sect. 3.2.2 of the electronic version of the paper). 
Then the scale height in the Solar neighborhood, scaled 
to r Q = 8.5 kpc, is hth = 600 pc; this value was taken for 
6-9 kpc in M51. 

For the outer regions of M51 we scaled the values de- 
rived from a comparison of Hn and Hi observations of 
the Milky Way. The Hi layer in the Milky Way becomes 
thicker at large radii, and the same may hold for the ion- 



ized gas. Dickey & Lockman (1990) derived a constant Hi 
scale height of 165 pc (the cool plus warm components) be- 
tween 4 and 8 kpc, and it increases considerably beyond 
the Solar circle (Henderson et al. 1982). Near the Sun the 
Hi scale height is about 200 pc, a factor of 3 smaller than 
the scale height of the ionized gas. Assuming this ratio to 
be constant for r > r Q , and taking the solar-neighborhood 
values at a radius of 9 kpc in M51, we found the scale 
heights of the thermal gas for the outer rings as given in 
Table 2. A sketch of M51 illustrating the spatial distribu- 
tion of the various components is shown in Fig. 2. 

The filling factor of the free electrons in the thermal 
disk of M51 was also adopted from the Solar neighbor- 
hood. It was calculated by comparing the flux density 
of the thermal radio emission with that expected from 
the disk of thermal electrons (Berkhuijsen, in prepara- 
tion). The thermal flux density at 2.8 cm of the Solar neigh- 
borhood, scaled to the distance of M51 and seen with 
a beam width of 75", is S2.8 = 0.35 ± 0.10 mJy/beam. 
Using Eq. (1), with C taken for the distance of M51, 
(n e ) = 0.035 ± 0.005 cm -3 and h th = 600 pc (as appli- 
cable near the Sun) we obtained / = 0.075 ± 0.030. This 
value is in good agreement with other estimates, being 
halfway between the filling factor of the diffuse warm gas 
(/ > 0.2, Reynolds 1991b) and that of giant Hn regions 
(/ ~ 0.01, Giisten & Mezger 1983). 

As no information is available on how / varies with r 
in the Milky Way, we adopted / = 0.075 for all 4 rings in 
M51. 
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Fig. 2. Sketch of the various layers in 
M51 seen in a vertical plane through the 
centre: disk of thermal electrons with 
scale height h t ^, synchrotron disk at 
A20.5cm with scale height and 
the elliptical halo. The depth in the 
thermal disk Az from which polarized 
emission at A20.5 cm is observed, is in- 
dicated. Z is the extent of the halo at 
z > hth- For clarity the thermal and 
synchrotron disk are shown thicker than 
to scale. 



3.3. Depolarization 

In Sect. 3.3.1 and 3.3.2 of the electronic version we dis- 
cuss wavelength-independent depolarization and Faraday 
depolarization outside the disk of M51. Based on that 
discussion, we assume that the Faraday depolarization 
occurs entirely in the thermal disk of M51. In order to 
take out the effect of the wavelength-independent depo- 
larization we shall use the relative depolarization between 
A = 20.5 cm and A = 2.8 cm, denoted as DP20.5/2.8 = 
P20.5fP2.8- As Faraday effects at A2.8cm are negligible, 
this ratio is essentially a measure of the Faraday depolar- 
ization at the longer wavelengths. The observed values are 
given in Tabic 3. 

Table 3. The observed relative depolarization between 
A20.5cm and A2.8cm, DP 2 o. 5/2.81 the vertical extent Az 
of the upper layer of the thermal disk visible in polarized 
emission at AA18.0/20.5 cm, and the corresponding values 
of £ (D) and £< H ) at A = 20.5 cm calculated from Eq. (7). 
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The Faraday depolarization in the disk is caused by 
both differential Faraday rotation and internal Faraday 
dispersion. Using results of Burn (1966) with the values 
in Tables 1, 2 and 3 we find that in the disk of M51 each 
of these effects is strong enough to significantly depolarize 
the emission at A20/18cm. Due to internal Faraday dis- 



persion, only an upper layer of the disk is visible. As we 
estimate below, this layer is only about 200-300 pc deep 
at r = 3-9 kpc. It can be easily seen that depolarization 
due to differential Faraday rotation across this depth is 
relatively weak. 

The fact that only polarized emission from an upper 
layer is observed is evident from the much smaller rotation 
measures observed between AA20.5 and 18.0 cm (Horellou 
1990) than between AA6.3 and 2.8cm (Neininger 1992b), 
at which wavelenghts Faraday rotation is negligible. For 
the two inner rings RM(20.5/18.0) ~ 0.25 x RM(6.3/2.8). 
As the disk is transparent at short wavelengths, this in- 
dicates directly that only a part of it is seen in polarized 
emission at A w 20 cm. The visible depth in the disk is then 
estimated as Az ~ 0.25 x 2h th = 200-300 pc for r = 3- 
9 kpc. In the case of field reversals in the part of the disk 
invisible at AA20.5, 18.5 cm this value is an upper limit 
to Az. Horellou et al. (1992), using different arguments, 
also concluded that at AA20.5, 18.0 only the upper part of 
the polarized disk is observed, and Beck (1991) found the 
same for NGC 6946. 

As the data at AA6.3, 2.8 cm are not complete for the 
two outer rings, we cannot make the above comparison for 
r = 9-15 kpc. Instead we propose the following estimate 
of the minimum visible depth. Let us define Az as the 
depth in the thermal disk from which polarized emission 
is observed (see Fig. 2). Then the layer in the synchrotron 
disk, which produces the observed polarized emission at 
A20.5 cm, has the thickness (/120.5 — ^th) + Az. If no Fara- 
day depolarization occurred in the visible layer, then the 
fraction DP20.5/2.8 °f the polarized emission at A20.0cm 
would come from a layer DP20.5/2.8 x 2/120.5 deep. Since 
this depth must be equal to the former value, this yields 

Az = /l20.5(2DP 20 .5/2.8-l) + ^th • (4) 
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As some depolarization actually occurs within Az, the 
true visible depth must be larger than this. The values 
of Az thus calculated are given in Table 3. These values 
are remarkably close to the upper limits derived above 
from the RMs observed in the two wavelength ranges. 
We note that in the radial range 12-15 kpc the thermal 
disk is completely transparent to polarized emission at 
AA18.0/20.5cm. 

3.3.1. Qualitative analysis of Faraday rotation in a 
double-layer system 

The polarization angle of the polarized emission is given 

by 

ip = RM fg A 2 + RM A 2 + Vo , (5) 

where RMf g is the foreground Faraday rotation measure 
produced mainly within the Milky Way. RM is the intrin- 
sic Faraday rotation measure produced by the magnetic 
field within the galaxy considered; A is the wavelength 
and ipQ is the intrinsic polarization angle. 

In order to distinguish between the contributions of 
the disk and the halo (see Sect. 4.1), we write 

RM = C (D) RM (D) + C (H) RM (H) , (6) 

where RM (D) and RM (H) are the Faraday rotation mea- 
sures produced across the disk and the halo, respec- 
tively, if both arc fully transparent to polarized emission. 
RM (D) = 0.81(n c )( D )B^ D) /i t h is defined here in terms of 
the disk scale height. (The full thickness of the disk is 
2/ith-) However, it is more convenient to define RM^ H ' = 
0.81(n o )( H )£^ H 'z in terms of the vertical extent of the 
halo, Z, i.e., the distance along z between z ~ hth and the 
upper boundary of the halo (see Fig. 2). B» is the line-of- 
sight component of the regular magnetic field. Since not 
the whole disk (or even halo) may be visible in polarized 
emission at a given wavelength, we introduce factors ^ D - ) 
and £( H ). As follows from above, £( D ) and £( H ) depend 
on the wavelength. We assume that both £( D ) and 
are the same at A2.8 and 6.2 cm, and also at A18.0 and 
20.5 cm. 

We can see from Tables 1 and 2 that the synchrotron 
disk at AA2.8 and 6.2 cm is about as thick as or thinner 
than the thermal one. Thus there is only little synchrotron 
emission originating in the halo and the halo magnetic 
field can be detected mainly via Faraday rotation in the 
near half. At AA18.0 and 20.5 cm, where h syn > h t h, the 
disk is not transparent to polarized emission at r < 10 kpc 
where the halo is present. As a result, it is impossible to 
determine the structure of the magnetic field in the part of 
the halo lying beyond the thermal disk from observations 
of the intrinsic polarized emission. 

Now we express and in terms of the scale 
heights of the thermal and synchrotron disk, h t h and h syn , 
and Az in a given wavelength range. One should take 



into account that, if synchrotron emission and Faraday 
rotation occur in the same region, the observed Fara- 
day rotation measure of a transparent layer is equal to 
\ J B\\n c dL, whereas that produced in a foreground Fara- 
day screen (i.e., a magneto- ionic layer devoid of relativistic 
electrons) is / B\\n c dL. Assume that h syn > hth, which 
inequality is true in the case of M51 for A w 20 cm. 
The Faraday rotation measure observed from the disk is 
given by iRM ( - D - ) Az / hth- The contribution of the halo to 
the observed Faraday rotation measure is iRM< H )(/i syn - 
h th )/Z + RM (H) [Z - (h syn - hth)]/Z, where the first term 
is due to the synchrotron-emitting region, and the second 
one is the contribution of the rest of the halo which acts 
as a foreground Faraday screen. Thus, we have 

£(D) = lA£ ; g (H) =1 _l fes y n-feth forZ>0 . (7) 

hth Z 

£W is undefined for Z = whereas RM (H) = in this 
case. Expressions for Az applicable in other cases can be 
found in Sect. 3.3.3 of the electronic version of the paper. 

At A = 2.8/6.2 cm the galaxy is transparent, so that 
Az = 2/i t h and we obtain £( D ) = 1; furthermore, £^ H ) s=s 1 
because at these wavelengths h syn differs insignificantly 
from hth for r = 3-9 kpc and h syn < h t h for r = 9 
12 kpc. However, £( D ) strongly differs from unity at A = 
18.0/20.5 cm. The dependence of £ (D) and on A is due 
to the A-dependence of Az and h syn . The values of Az, 
£( D ) and given in Table 3 refer to AA18.0/20.5 cm. 
The halo is transparent for polarized emission at all the 
wavelengths considered, and £( H ) differs from unity only 
because some synchrotron emission originates within the 
halo (at h syn < z < hth), whereas the remaining part of 
the halo acts as a foreground screen. 

4. Recognition of magnetic field patterns 

When fitting the observed distribution of the polariza- 
tion angle ip at a given wavelength A we adopt for each 
layer the following truncated Fourier representation for 
the cylindrical components of the regular magnetic field 
-E?rcg = (B r , Bg,B z ): 

B r = B sinpo + Bi sinpi cos(# — (3) , 

Bg = B cospo + B\ cospi cos(6* — (i) , (8) 

B z = B za + B zl cos(9~ (3 Z ) , 

where B and B z0 are the strengths of the horizontal (par- 
allel to the galactic plane) and vertical (perpendicular to 
the plane) components of the m = mode, respectively, 
Bi and B zi are those of the m = 1 mode, p and p\ are 
the pitch angles, and (3 and (3 Z are the azimuthal angles 
at which the corresponding non-axisymmetric components 
are maximum. In Eq. (8) only the two lowest modes have 
been retained; this proves to be sufficient to fit the avail- 
able data. The magnetic pitch angle is the small angle 
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measured from the magnetic field vector to the tangent 
of the local circumference. It is positive (negative) if the 
magnetic field spiral opens counterclockwise (clockwise). 
We note that the magnetic field direction can be either 
inwards or outwards along the spiral. In the case of M51 
a negative pitch angle corresponds to a trailing spiral. 

The intrinsic polarization angle tpo in Eq. (5) is de- 
termined by the transverse component of the magnetic 
field. A suitable expression relating Vo to the magnetic 
field of the form given in Eq. (8) was given by Sokoloff 
et al. (1992); its detailed derivation can be found in Ap- 
pendix A of the electronic version - see Eq. (A3). As the 
synchrotron emissivity in the halo is significantly weaker 
than in the disk, we assume that the intrinsic polarization 
angle depends solely on the field in the disk. 

We should emphasize that we analyze simultaneously 
and consistently the longitudinal and transverse (with re- 
spect to the line of sight) components of the magnetic field 
which manifest themselves through RM and tpo, respec- 
tively. Previous work either considered solely the Faraday 
rotation measures between pairs of wavelengths or simpli- 
fied the model by supposing that Vo = const (see Sokoloff 
et al. 1992). In both cases only the line-of-sight magnetic 
field could be recovered from observations. Attempts to 
extract additional information about the transverse com- 
ponent of the magnetic field from an independent analysis 
of "magnetic pitch angles," or ipo, from total and polar- 
ized intensity, etc. often led to results that were inconsis- 
tent with those obtained from the RM analysis. Here for 
the first time we propose a consistent three-dimensional 
model of the regular magnetic field observed in a galaxy. 

The fitting procedure and the estimation of uncertain- 
ties of the fitted parameters are discussed in detail in Sect. 
4.2, 4.4, and Appendix B of the electronic version of the 
paper. 

In order to obtain satisfactory fits to the data, we cal- 
culated the residual S defined in Eq. (Bl), that charac- 
terizes the deviation of the fit from the measured points, 
found its minimum with respect to the above fit parame- 
ters, and then employed the x 2 an d Fisher statistical tests 
briefly discussed in Appendix B to assess the reliability of 
the fit. The \ 2 statistical test, Eq. (B2), ensures that the 
fit is close enough to the measured points with allowance 
for their weights equal to <r~ f with a n i the standard devi- 
ation of ipni (see Sect. 2). The Fisher test, Eq. (B3), was 
then applied to verify that the quality of the fit is the same 
at all individual wavelengths (see Sokoloff et al. 1992 for 
details). 

4-1. Statistical evidence for a magneto-ionic halo in M51 

Fits for all wavelengths obtained with one transparent 
magneto-ionic layer were generally inconsistent with the 
Fisher test at r < 9kpc. Therefore, we tried fits to 
the polarization angle distributions for AA2.8/6.2 cm and 
AA18.0/20.5 cm separately. The magnetic patterns re- 



vealed in the two wavelength ranges in the inner two 
rings turned out to be very different from each other. 
For example, the values of B obtained at the shorter 
and longer wavelengths had different signs. Apparently the 
only plausible explanation for this is that different regions 
are sampled in the two wavelength ranges and that the 
magnetic field has completely different configurations in 
them. Moreover, Bu changes sign either in the upper part 
of the disk Az or in a layer above z = h t h- Since it is very 
unlikely that the regular magnetic field may have such a 
complicated vertical structure within the upper part of 
the disk, we concluded that the change of direction of 
suggests the existence of one more extended component in 
the galaxy. Thus, at r < 9kpc, at least two extended com- 
ponents of the magneto- ionic medium in M51 are present, 
namely the disk and the halo. Then the observed Fara- 
day rotation at the short wavelengths is dominated by the 
magnetic field in the disk whereas that at the long wave- 
lengths is mainly determined by the field in the halo. 

It is clear that separate fits for different wavelength 
ranges cannot be physically satisfactory. Therefore all the 
fits discussed below were obtained from all four wave- 
lengths simultaneously using the values of Az in Table 3. 
For the inner rings, r — 3-9 kpc, we applied a double-layer 
model of the magneotionic medium using the values of Az 
in Table 3, whereas for the radial range 9-15 kpc the fits 
were performed for a single-layer model without distin- 
guishing the disk and halo contributions to RM (that is, 
we put B\ =0 there). 

Evidence of a halo in M51 from X-ray observations is 
discussed in Sect. 4.3 of the electronic version of the paper. 
We adopted an elliptical shape for the halo with a height 
of 6 kpc above the midplane near the center (see Fig. 2) 
and assumed an electron density of 3 x 10 ~ 3 cm~ 3 as in 
the Milky Way. 

4-2. Results of the fitting 

The results of the fitting are presented in detail in Sect. 4.4 
of the electronic version of the paper. Here we combine 
parts of Tables 4 and 5 and we only mention that the 
regular magnetic field in the disk is shown to be a super- 
position of axisymmetric and bisymmetric modes with a 
slight predominance of the latter. The field in the halo ap- 
pears to be axisymmetric. An example of the fits is shown 
in Fig. 3 for r = 3-6 kpc. Similar figures for the other rings 
can be found in the electronic version. 

5. Global structure of the regular magnetic field 

In this section we derive the strength and the direction 
of the regular magnetic field, and we discuss the global 
properties of these parameters. 

The results are compiled in Tables 4 and 5. A note 
of caution is appropriate here: the resulting amplitudes 
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Tables 4 and 5 (combined). A fitted model for the regular magnetic field in M51 



r [kpc] 3-6 6-9 9-12 12-15 
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Disk 
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[PG] 
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6 9+ 3 ' 4 


3.4 ±o.8 




i o+ 2 - 6 


1 Q+9-0 
± - a -1.3 


4 4+40.2 

-3.2 



Notes: 

' a ' Dots mean that the corresponding parameter is insignificant. Errors are 2a Gaussian deviations. 
^ RMfg for 12-15 kpc was fixed at the given value (see Sect. 4.4.4 of the electronic version). 



Fig. 3. The polarization angles (dots with error bars, measured from the local radial direction in the plane of M51) as a function 
of azimuthal angle in the galactic plane and fits (solid lines) for the ring 3 < r < 6 kpc at (a) A2.8, (b) A6.2, (c) A18.0 and (d) 
A20.5 cm. Error bars show the la errors of the measurements. Dashed lines indicate the foreground levels, RMf g A 2 
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B r eg(n e }Z for the halo with scaling factors specified in 
the caption. 

In the outer ring, 12-15 kpc, the magnetic field struc- 
ture is distorted. The values of p and p\ differ consider- 
ably from those for r < 9 kpc and are even positive. In- 
spection of the polarization map in Fig. la confirms that 
in the northern part the magnetic pattern at these radii 
is plagued by strong distortions still having a rather large 
spatial scale. Only in the ring 12-15 kpc the magnetic field 
has a weak vertical component. 



Fig. 7. The radial variation of the strength of the regular mag- 
netic field obtained from our fits and averaged in the rings (cir- 
cles with error bars), the total magnetic field obtained from 
the total intensity of the nonthermal emission assuming en- 
ergy equipartition between cosmic-ray particles and magnetic 
field (dashed) and its regular component obtained using the 
observed degree of polarization (solid). 

Bi were obtained assuming that (n e ) is independent of 
azimuthal angle. For each ring we also give the strength 
of the regular magnetic field averaged over the azimuth, 
B lcg . Uncertainties were calculated with allowance for er- 
rors in (n c ) too. The large errors, especially at 12-15 kpc, 
are mainly due to the errors in (n c ). 

In Fig. 7 we show the radial variation of B leg and that 
of the total and regular magnetic field strengths, B and 
B rcg , obtained from the total nonthermal emission and the 
observed degree of polarization as described in Sect. 3.1 
and given in Table 1. We note that for each ring there 
is a close agreement between the two values of the reg- 
ular magnetic field which were obtained from completely 
independent physical parameters and methods. Between 
r = 3 and 15 kpc the exponential radial scale length is 
14.3 ±6.0 kpc. 

As can be seen from Tables 4 and 5 the magnetic fields 
in the halo and in the disk inside r = 12 kpc are hor- 
izontal. In the halo we have Bq ' > 0, p y ' < and 
B\ n > ~ 0, which means that the radial component of 
the regular magnetic field is directed inwards, with the 
azimuthal component directed counterclockwise; that is, 
Bl H) < and B { e H) > - see Eq. (8). Meanwhile, in the 
disk we have ^ D) w b[ d) and both are negative together 
with the pitch angles. Therefore, for r = 3-9 kpc the ra- 
dial field in the disk is directed outwards, with Bg directed 
clockwise at almost all 9. We conclude that the regular 
magnetic fields in the disk and the halo have almost oppo- 
site directions everywhere within r = 9 kpc except in the 
northwestern part of the ring at r = 3-6 kpc. 

In Fig. 8 the direction of B lcg in each sector is shown 
for the disk and the halo separately. The length of the 
vectors is proportional to B Icg (n c )h t h for the disk and 



6. Discussion 

6.1. The magnetic field in the halo 

The available polarization measurements performed at the 
two pairs of widely separated wavelengths allowed us to 
determine the magnetic field structure in two regions along 
the line of sight. In the text above we called these regions 
the disk and the halo. This usage was justified in Sect. 4.1. 

Our results represent the first indication of a magneto- 
ionic halo in a galaxy seen nearly face-on. The detection 
of a radio halo in an edge-on galaxy is a difficult obser- 
vational problem, even more so the determination of the 
magnetic field structure. In the galaxies seen nearly face- 
on some of the difficulties are alleviated. First, the halo 
and its magnetic field are illuminated by a strong back- 
ground source of polarized emission, the disk. Second, the 
polarization measurements over the entire disk can be used 
to reveal the global azimuthal structure of the field in the 
halo as it was done in the present paper. It is important 
to note that the magnetic field in M51 has different struc- 
tures in the disk and the halo. If the field structures were 
similar to each other, the detection of the magneto-ionic 
halo might be difficult. 

We showed that the field in the halo of M51 is pre- 
dominantly horizontal as in the halos of NGC 891 and 
NGC 253 (Hummel et al. 1991b; Sukumar & Allen 1991; 
Beck ct al. 1994). 

According to our fits we estimate the halo radius to be 
about 10 kpc. This estimate agrees with the data on X- 
ray emission from M51 which also indicate a halo radius 
of about 10 kpc (Ehle ct al. 1995). 

With the values of (n e ) and Z from Table 2, the esti- 
mated strength of the regular magnetic field in the halo 
decreases from about 3 /zG at the radial distance of 3- 
6 kpc to about 1 fiG at r — 6-9 kpc. The field is basically 
axisymmetric. The upper limits on the m = 1 mode in 
the halo arc estimated from our fits as \B[ n, \ < 1/iG and 
< 3 /iG for r = 3-6 and 6-9 kpc, respectively. 

It is interesting to compare the values of the regular 
magnetic field in the halo with the upper limit on the total 
magnetic field strength estimated from the equilibrium be- 
tween thermal and magnetic energy densities in the X-ray 
emitting gas (Ehle et al. 1995). With (n c ) = 0.003 cm- 3 
and a volume filling factor of 0.8, their results yield 
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Fig. 8. Directions of the horizontal regular magnetic field in the disk (a) and halo (b) of M51 according to the fits presented in 
Tables 4 and 5. For clarity we scaled the vectors as follows: they are proportional to 1Z (= B Icg {n c }h t h) in the inner two rings 
in the disk, to 37?. in the two outer rings in the disk, to 2.571^ (= B les {n e )Z) for 3-6 kpc in the halo, and to 571^ for 6-9 kpc 
in the halo. The vertical component at r = 12-15 kpc was not included. The vectors are shown superimposed onto an optical 
picture of M51. The sectors and rings used are indicated. 



B < 7 /iG. Our results are consistent with this limit. If 
the true total field strength is close to the above upper 
limit, the turbulent field in the halo has a strength of 
about 6 /uG exceeding that of the regular magnetic field. 

The global field directions are in general opposite in 
the disk and the halo of M51. This implies that the regu- 
lar magnetic field in the halo cannot be simply advected 
from the disk. Such reversals appear in the dynamo the- 
ory for galactic halos (Ruzmaikin et al. 1988, Sect. VIII. 1; 
Sokoloff & Shukurov 1990; Brandenburg et al. 1992) and 
could be due to the topological pumping of magnetic field 
by a galactic fountain flow (Brandenburg et al. 1995). 
Moreover, the dominance of the axisymmetric field in the 
halo is also consistent with the mean-field dynamo theory 
which predicts that non-axisymmetric magnetic modes 
can be maintained only in a thin galactic disk and most 
likely decay in a quasi-spherical halo (see Ruzmaikin et 
al. 1988). We cannot say anything about the parity of the 
halo field with respect to the midplane because the galaxy 
is not transparent at A « 20 cm in the rings where the halo 
is present and, in addition, the synchrotron emission from 
the halo is negligible. 

6.2. The azimuthal structure of the field 

The azimuthal distributions of polarization angle in M51 
seen over the radial range r = 3-15 kpc are success- 



fully represented by a superposition of only two azimuthal 
modes of the large-scale magnetic field, m = and m = 1 
in the disk. Even though we restrain ourselves from iden- 
tifying these magnetic harmonics with dynamo-generated 
axisymmetric and bisymmetric modes before a more care- 
ful theoretical analysis has been made, we mention that 
dynamo theory also predicts that the two leading az- 
imuthal modes m = and m = 1 typically dominate in 
spiral galaxies. Furthermore, it follows from the dynamo 
theory that non-axisymmetric magnetic structures should 
be more often a superposition of the two azimuthal modes 
than a purely bisymmetric mode (Ruzmaikin et al. 1988, 
p. 231; Beck et al. 1996). A similar superposition of modes, 
but with a dominance of the bisymmetric mode, was found 
earlier in M81 by Sokoloff et al. (1992) (see also Krause 
et al. 1989b). In M31 the axisymmetric magnetic mode is 
dominant (Ruzmaikin et al. 1990), and for NGC 6946 a 
superposition of m = and m = 2 magnetic modes was 
suggested by Beck & Hoernes (1996). 

Of course our results do not imply that higher az- 
imuthal magnetic modes are not present in M51, but only 
that the accuracy of the available observations is insuffi- 
cient to reveal them. One can expect that the amplitudes 
of the harmonics with m > 2 arc considerably smaller 
than those of the modes m — and 1. 

Since the theory of the galactic mean field dynamo 
predicts an efficient generation of the bisymmetric mode 
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in M51 with the maximum of the m = 1 eigenmodc at r ~ 
2kpc (Baryshnikova et al. 1987; Krasheninnikova et al. 
1989), we are tempted to identify the m — 1 mode revealed 
for 3 < r < 9kpc with a bisymmctric field generated by 
the dynamo. This suggestion is confirmed by the closeness 
of the pitch angles p and p\ of the m = and 1 modes to 
each other in the two innermost rings: this is typical of the 
dynamo-generated fields in a thin disk (Ruzmaikin et al. 
1988). This conclusion is also plausible for r = 9-12 kpc. 
We also note that a nonlinear dynamo model of Bykov et 
al. (1996) predicts a mixture of magnetic modes, which 
is roughly similar to that detected here, to be found in 
some vicinity of the corotation radius, i.e. just near 6 kpc 
in M51. 

The azimuthal modes inferred for the outermost ring 
can be hardly identified directly with the dynamo modes 
because the pitch angles of individual modes are posi- 
tive. These modes may be due to distortions imposed by 
non-axisymmetric density and velocity distributions pos- 
sibly caused by the encounter with the companion galaxy 
NGC 5195 (Howard & Byrd 1990). Concerning the total 
horizontal regular magnetic field, its pitch angle is nega- 
tive for 70° < 9 < 170° and positive in the rest of the 
sectors. Inspection of polarization maps in Fig. 1 confirms 
that the pattern of polarization angles at these radii in 
the northern part is strongly distorted on a rather large 
scale. 

We note that the regular magnetic field in the disk of 
M51 is directed outwards, whereas those in IC 342 (Krause 
et al. 1989a), M31 and NGC 6946 (Beck et al. 1996) are 
directed inwards. 

6.3. Inner and outer spiral structure 

Inspection of Tables 4 and 5 shows that for the disk most of 
the fitted parameters of the inner rings (r < 9 kpc) differ 
systematically from those of the outer rings. The phase 
angle [3 varies by about 110° between the rings at 3-6 kpc 
and 6-9 kpc, and also between 9-12 kpc and 12-15 kpc, but 
not between 6-9 and 9-12 kpc. The inner pattern is more 
coherent and has stronger magnetic field than the outer 
pattern. Thus it seems that the magnetic field structure 
in the outer rings is not a smooth continuation of the 
structure in the inner rings, but that rather two distinctly 
different magnetic field structures are present in M51. 

This result is very interesting as Elmegreen et al. 
(1989) showed, using optical plates, that M51 contains 
an inner and an outer spiral structure which are overlap- 
ping between r = 6 kpc and r = 8 kpc (see Sect. 6.3 in 
the electronic version for a more detailed discussion) . The 
outer spiral arms are thought to be material arms driven 
by the companion, whereas the inner spiral arms are due 
to density waves caused by the outer arms. 

A discontinuity in the magnetic field pattern indicates 
that different physical effects contribute to the field struc- 
ture at r < 9 kpc and r > 9 kpc. The discontinuity occurs 



at the radius where the inner spiral structure ends and 
the outer spiral structure becomes dominant. Therefore, 
the two magnetic field structures can be physically con- 
nected with the inner and outer spiral patterns proposed 
by Elmegreen et al. (1989). 

The relatively strong magnetic field and its regular 
pattern in the inner region are compatible with the idea 
of a dynamo acting under more or less steady conditions. 
In the outer regions, where the spiral arms are produced 
by a recent encounter with NGC 5195 about 10 8 years 
ago (Howard & Byrd 1990), the magnetic field may be 
a remnant of an older one disrupted by the velocity per- 
turbation. Therefore it is understandable that the pitch 
angles are irregular and the magnetic field is weak. 

We conclude that the magnetic field pattern in the disk 
of M51 appears to be not one global structure, but consists 
of an inner pattern associated with the inner spiral struc- 
ture of density wave arms and an outer pattern related to 
the outer spiral structure of material arms. The interac- 
tion between magnetic fields and the spiral patterns is not 
yet understood. 

6.4- Pitch angles of the magnetic field and of the spiral 
arms 

We compared the pitch angles of the magnetic field in the 
disk derived from Eq. (10) (see Sect. 6.4 in the electronic 
version) with the pitch angles of the dust lanes running 
along the inner edge of the optical spiral arms as tabulated 
by Howard & Byrd (1990). In each ring the optical pitch 
angles were averaged in the same sectors as were used for 
the model fits (see Fig. 8). The comparison was possible 
only for the inner two rings, as at larger radii the measured 
optical pitch angles and the magnetic model pitch angles 
have too few sectors in common. 

Comparing the corresponding sectors we found general 
agreement between optical and magnetic model pitch an- 
gles. For the ring 3-6 kpc the mean of the optical pitch 
angles is —15° ± 8° and that of the magnetic pitch angles 
is —11° ±3°, whereas for the ring 6-9 kpc these values are 
-13° ± 12° and -10° ± 8°, respectively. The errors are 
one standard deviation from the mean value and are due 
to intrinsic variation in pitch angle in each ring. Although 
the agreement is quite good, we note that the optical pitch 
angles show larger variations than those of the fitted mag- 
netic field. 

Altogether, we conclude that on average the magnetic 
field inferred from our fits is well aligned with the spiral 
arms, although local misalignments may be considerable 
(see Fig. 8). 

Section 6.5 The origin of the vertical field and also Ap- 
pendices A. The intrinsic polarization angle, B. Statistical 
tests and errors and C Basic notation can be found in the 
electronic version of the paper. 
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7. Conclusions 

1. The global magnetic pattern in M51 at 3 < r < 15kpc 
can be represented as a superposition of the two lowest 
azimuthal Fourier modes. 

2. Our analysis indicates the existence of a magneto- 
ionic halo and shows that the magnetic fields in the disk 
and the halo have different configurations. 

3. The radial extent of the halo is about fOkpc, in 
agreement with X-ray data. The halo field is horizontal 
and axisymmetric. The regular magnetic fields in the halo 
and in the disk are spirals. The field directions along the 
spirals are generally opposite running inwards and out- 
wards in the halo and the disk, respectively. 

4. In the disk a superposition of axisymmetric (m = 
0) and bisymmetric (m = 1) magnetic modes provides 
a satisfactory fit to the observations. The m = 1 mode 
slightly dominates at 3 < r < 6kpc and the two modes 
have about equal amplitudes at 6 < r < I5kpc. The field 
is predominantly horizontal between 3 and 12 kpc and has 
a weak vertical component at 12 < r < 15 kpc. In the 
rings between 3 and 9 kpc the field structure in the disk is 
strongly non-axisymmetric with a field maximum in the 
eastern part of the galaxy and a weak field in the western 
part. Details are given in Tables 4 and 5 and shown in 
Fig. 8. 

5. The magnetic field pattern in the disk of M51 shows 
a discontinuity at r ~ 9 kpc, the position at which the 
inner and the outer spiral structure join (Elmcgrcen ct 
al. 1989). The relatively strong, coherent magnetic field in 
the inner rings occurs in the system of spiral arms excited 
by density waves, whereas the weaker and partly distorted 
field in the outer rings exists in the area of the material 
spiral arms produced by the encounter with the compan- 
ion. 

6. The azimuthally averaged strength of the regular 
magnetic field obtained for the disk decreases from about 
7/xG at radius 3-6 kpc to about 4,uG at 12-15 kpc (see 
Fig. 7). The strength of the regular field in the halo de- 
creases from 3 [iG at 3-6 kpc to zero beyond 9 kpc. Details 
are given in Tables 4 and 5. 

7. The azimuthal averages of the regular magnetic field 
strength in the disk obtained from our fits are in good 
agreement with independent estimates from the total syn- 
chrotron emission and the degree of polarization. The ra- 
dial scale length of the regular magnetic field is 14 ±6 kpc. 

8. We compared the pitch angles of the regular mag- 
netic field obtained from the fits with those of the dust 
lanes delineating spiral arms for the sectors they have in 
common. Their mean values agree to within the errors. 
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